Phase Estimation With Interfering Bose-Condensed Atomic Clouds 
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We investigate how to estimate from atom-position measurements the relative phase of two Bose- 
Einstein condensates released from a double-well potential. We demonstrate that the phase es- 
timation sensitivity via the fit of the average density to the interference pattern is fundamentally 
bounded by shot noise. This bound can be overcome by estimating the phase from the measurement 
of yN (or higher) correlation function. The optimal estimation strategy requires the measurement 
of the iV-th order correlation function. We also demonstrate that a second estimation method - 
based on the detection of the center of mass of the interference pattern - provides sub shot-noise 
sensitivity. Yet, the implementation of both protocols might be experimentally challenging. 
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Introduction. The experimental realization of the 
Bose-Einstein condensation (BEC) opened a new chapter 
in the field of atom interferometry [l| . The BEC consti- 
tutes a bright and well-controllable source of particles, 
in analogy to the laser light commonly used for interfer- 
ometric purposes. Atoms are very good candidates for 
precise measurement of electromagnetic [2-4] or gravita- 
tional [a, [6j forces. Moreover, non-classical states have 
been created via Quantum Nondemolition measurements 
[7|, |8| , or exploiting the naturally present two-body inter- 
actions [9l4llj. These states allow to overcome the limit 
imposed by the classical physics on measurement preci- 
sion (shot- noise limit) [12|, [l3|, as recently demonstrated 
experimentally }l4J . An atom interferometer can be im- 
plemented using a gas trapped in a double-well poten- 
tial [15l4l7| , where the two spatial modes localized in its 
minima play the role of the interferometer's arms. In the 
simplest interferometric sequence, first a relative phase 
9 builds up between the wells and then the BECs are 
released from the trap and form an interference pattern. 
In this manuscript we analyze different phase estimation 
strategies relying on position measurement of the atoms 
forming this pattern. We show that fitting the average 
density to the interference pattern [15j , gives a sensitivity 
A9 bounded by the shot-noise. This limit can be over- 
come by measuring correlations between the atoms of 
order not smaller than y/N , where N is the total number 
of atoms. We also demonstrate that the estimation by 
the N-th order correlation function is optimal - the sen- 
sitivity saturates the bound set by the Quantum Fisher 
Information (QFI) 18]. Finally, we present another esti- 
mation scheme which can give sub shot-noise sensitivity, 
based on the measurement of the position of the center 
of mass of the cloud. Although it involves a probabil- 
ity which is simpler to construct with respect to high- 
order correlation functions, it requires detection of all N 
atoms. These results indicate that the achievement of 
sub shot-noise sensitivity using the interference pattern 
might prove challenging. 

The model. We introduce the two-mode field operator 
of a bosonic gas in a double- well potential, ^>(x,t) — 



ip a (x,t)a + ip b (x,t)b, where a) /W creates an atom in 
the left/right well. While the atoms remain trapped, 
an unknown relative phase 9 is imprinted between the 
wells. This stage is described by a unitary evolution 
U{9) = e~ l6J * of the initial state \ipi n ) of the double- 
well system [19|. After the phase acquisition, the trap 
is switched off and the initially localized mode functions 
ip a / b (x,t) freely expand. 

All information which can be extracted from a posi- 
tion measurement is contained in pn{xn\0) [2CH22|- the 
conditional probability of finding N particles at posi- 
tions xn — {xi . . . xn). It can be expressed in terms 
of the N-th order correlation function [23| as follows: 
Pn{xn\9) = tftGn(3n,0)- To evaluate this probabil- 
ity, we write the initial state in the well-population ba- 
sis, \ipin) — J2n=o C n \ n i N~ n) and suppose that the 
expansion coefficients are real and posses the symme- 
try C„ = Cjv-n 24]. It is convenient to use the 
Heisenberg representation and evolve the field operator, 
U^(9)^(x,t)U(9) = ii a (x,t)e l! la + ip b (x,t)e~^b. Then, 
we expand the Fock states \n, N — n) in the basis of the 
spin-coherent states [25| and obtain 
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where u e (x,ip;t) = ip a (x,t)ei^ +d ^ +ip b (x,t) e -i^ +e K In 
the remaining part of the manuscript, the time t is cho- 
sen such that the interference pattern is already formed. 
The probability ([1]) is the starting point for the following 
discussion of various phase estimation strategies. 

Estimation via the fit to the density. One of the most 
natural ways of estimating the value of 9 is through 
the fit of the average density to the interference pat- 
tern. The experiment consists in dividing the inter- 
ference pattern into M bins located at positions Xi 
and measuring the number of particles n% in each of 
them. This procedure is repeated to times: the set 
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n%"', k = l,.,.,m is recorded and the average rij 



X^feLi n i / m m *"th bin is evaluated. The average oc- 
cupation (rii) = lim m _ i . 00 X^feLi n i V 771 ( we assume to be 
known a priori in the experiment) with free parameter 9 
is then fitted to the histogram of the measured density 
{xi,fii},i — 1, ...,M. For small bin size Ace, (nt) is re- 
lated to the average density by (n,) = G\(xi,9)Ax. The 
value of 9 can be determined from the least squares for- 

0. The number fluctuations 
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mUla d§ 2^i=l 2A i n,/m 

A 2 ni = lim m _ i . 00 ^feLi^i"'' — ( n i)) 2 (also assumed to be 
known a priori in the experiment) can be calculated from 
the probability p(rii \ 9) of detecting rii particles in the «-th 
bin, 
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Using Eq.((T]) we obtain 



A 2 n t = G 1 (x l ,9)Ax + [G 2 (x t , Xl ,9) - G\{ Xi ,d)] (Ax) 2 . 

(2) 
In order to derive the sensitivity A 2 9 for this esti- 
mation method, we notice that, if the number of mea- 
surements to is large, according to the central limit 
theorem the probability distribution for the average 
rii in the z-th bin tends to the Gaussian p{fii\&) — 

1 
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— _, . Although in every shot the atom 

counts are correlated between the bins, averaging over 
large number of experiments washes out this dependence. 
Therefore, the total probability of measuring the values 
{n} = (ni . . . n M ) is a product P({n}\9) = Yifii Pi^ild) ■ 
For this Gaussian probability, the condition for the least 
square fit coincides with the condition for the maximum 
likelihood estimator (MLE), ^P({n}\6) = 0. This is a 
crucial observation since it is well known [26[ that, for 
this choice of the estimator, the sensitivity is given by 
the Cramer- Rao Lower Bound, A 2 9 = F _1 . Here, F is 
the Fisher information (FI), 



F 



N 

E 



i 



d 



P({n}\9) \d9 

n 1 ...fi M =Q VL Jl ; \ 



P({n}\9) 



m>l 



->■ TO 



M 1 



i=l 



d(nj 
A 2 n t V 89 



(3) 



simplifies to 
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with the one-particle probability pi(x\9) = 
l(\Mx,t)\ 2 MMx,t)\ 2 )^(J x )Re[r a (x,t)Mx,t)e ie ]. 
It can be demonstrated [27j that Eq.(U]) gives F < mN, 
thus A9 > A9sn = ~m~m f° r an V input state (A9$n 
denotes the shot-noise sensitivity). Below we argue that 
the inclusion of the second term in the fluctuations does 
not improve the sensitivity. 

In Eq.@, the first term Gi(xi,9)Ax scales linearly 
with TV while the second term, as a function of N, is 
a polynomial of the order not higher than two, aN 2 + 
bN + c. Since the fluctuations A 2 n^ must be positive, 
then a > 0. Otherwise, for large N, no matter how small 
Ax, we would have A 2 n^ < 0. If a > 0, the positive sec- 
ond term enlarges the fluctuations and thus worsens the 
sensitivity. When a = 0, the first and second terms in 
Eq.@ scale linearly with N, and for small Ax the second 
term can be neglected, thus we end up again with Eq.Q 
(a similar argument can be made in order to conclude 
that increasing the bin size Ax always worsens the sensi- 
tivity with respect to Eq.Q). This is the first important 
result of this manuscript: A9sn is a lower bound for the 
fit sensitivity. This is because the maximal value of the 
FI is expressed in terms of the single-particle probability 
and does not exploit the correlations between the atoms. 
Let us now demonstrate that the measurement of these 
correlations can indeed improve the phase sensitivity. 

Estimation via the correlation functions. The estima- 
tion method we discuss now relies upon determination 
of the phase 9 from the fc-th order correlation function 
Gk{xk\9). We can choose to deduce 9 using the MLE 
[28|: given a fc-tuple £& of detected atoms' positions, 
the estimated value of the phase is chosen by maximiz- 
ing Gk(£,k\9) with respect to 9. The variance A 2 # after 
m 3> 1 experiments is given by A 2 9 = FT,} , where 
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Since conditions for the MLE and the fit coincide, the 
sensitivity of the latter is given by the inverse of <j3j> as 
well. 

We now demonstrate that this sensitivity is bounded 
by the shot-noise. Let us assume for the moment that 
the second term in the Eq. ([2]) - which is proportional to 
(Ax) 2 - can be neglected. If this is the case, the particle 
number distribution is Poissonian and the FI from ([3)) 



with pfc(xfc|#) = jv! Gk(x*k,9). The coefficient ^ ac- 
counts for the number of independent drawings of k par- 
ticles from TV, i.e. (^/(^Ti)- We notice that by setting 
k = I, i.e. the estimator is a single-particle density, we 
recover the FI from Eq.([3]). Therefore, the measurement 
of positions of N particles independently is, in terms of 
sensitivity, equivalent to fitting the average density, and 
thus is limited by the shot-noise. 
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FIG. 1: (color online). The sensitivity ^frnAQ (black solid 
lines) for N = 8 calculated by numerical integration of Eq. © 
for various A: as a function of \ifiin) £ A. The two limits, 
^/mA9sN and yfmAOHL are denoted by the upper and lower 
dashed blue lines, respectively. The optimal sensitivity, given 
by the QFI, is drawn with red open circles. The inset magni- 
fies the vicinity of the spin-coherent state, showing that the 
sub shot-noise sensitivity is reached starting from fc m i n = 4. 



Let us calculate the FI for the case k = N. As the 
interference pattern is formed after a long expansion 
time, the mode functions can be written as i/j a /b{ x Tt) — 

e 1 ^ 1 ^ 1 ° 2 • V> {jpz) , where a = \/—, "0 is a Fourier trans- 
form of the initial wave-packets (common for ip a and tpb) 
and the separation of the wells is 2xo- The above mode 
functions give the probability ((T|), which is inserted into 
Eq. © • After integration over space we get 



N 



\k=N) = rn 



*E C « 



n=0 



N 
~2 



m ■ 4AV 2 



(6) 



This is the second important result of this manuscript: 
the estimation of 9 from the iV-th order correlation func- 
tion is optimal - it saturates the QFI, which is a result 
of maximization of the FI with respect to all possible 
measurements [18j, [29j. According to Eq.©, all states 
for which A 2 J Z > -j give Ad < A9$n- Looking for sub 
shot-noise sensitivity, we take a sub-set of these states - a 
family of ground states of a double- well system with neg- 
ative interactions, as this is a natural choice in the con- 
text of this work. These states range from spin-coherent 
to the NOON state. For all states in this family, which 
we denote by A, the estimation by the measurement of 
the iV-th order correlation function gives A9 < A9$n, 
and reaches the Heisenberg limit A9ul = -i=-h for the 
NOON state (as indicated by open circles in Fig. [T]). 

It is important to know whether it is still possible to 
achieve sub shot-noise sensitivity by using lower order 
correlation functions k < N. As the space integrals in 
Eq.© for k 7^ N cannot be evaluated analytically, we 
calculate F numerically taking Gaussian wave-packets, 



with the initial width of the 
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wave-packets <7o = 0.1 and half of the well separation 
xq = 1. In Fig. Q]we plot the sensitivity using Eq.© for 
N = 8 atoms as a function of \ipi n ) £ A. We observe that 
the sensitivity improves with growing k and overcomes 
the shot-noise limit at fc m ; n = 4. 

Moreover, we numerically checked that, for larger N, 



kmin tends to v N (more precisely to the closest integer). 
However, measuring such a (typically) high order corre- 
lation function, during the required calibration stage ( pk 
must be known before performing the phase estimation 
protocol described above), involves a large experimental 
effort. Indeed, for every 9, the /c m j n -dimensional function 
j9fc min (xfc min , 0) must be probed. In the following section, 
we present a different detection scheme which gives sub 
shot-noise sensitivity and does not require the knowledge 
of a multi-dimensional function. 

Estimation via the center of mass measurement. Esti- 
mation of 9 from the measurement of the center of mass 
requires a relatively simple calibration stage. Positions of 
N atoms are recorded and from this data location of the 
center of mass is calculated. Such observable is described 
by the one-dimensional probability p C m(x\9), related to 
the full A^-body probability ©) by 



Pc,n(x\9) = / dx N S Ix - — ^ Xj\ 



p N (x N \9), 



where "S" is the Dirac delta. Modeling the mode- 
functions by Gaussians, as in previous section, and under 
the realistic assumption e~ x o/' J o <1, we obtain 



Pcm(x\9) (X e ~^ 
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Notice that the above probability depends on 9 only for 
states with non-negligible NOON component, i.e. non- 
zero Co (and thus Cjv) 30]. 

When p cm (x\9) is known, the phase can be estimated 
using the MLE analysis described in the previous section. 
The sensitivity is given by the inverse of the FI which 

can be calculated analytically using the definition F cm — 

00 

m / dx PaJ(x\6) (wPcm{x\9)) . We obtain 



F nm = mN z 



i .- a/1-2 Cl 



(7) 



In Fig [2] we plot the sensitivity calculated by the inverse 
of the FI ([7]) as a function of \ipi n ) € A. The estimation 
through the center of mass is not optimal, but the sen- 
sitivity can still be better than shot-noise, and tends to 
A9 HL for \ip in ) -> NOON. We underline that, although 
the calibration stage is not as difficult as in the case of 
high-order correlations, phase estimation based on the 
center of mass measurement requires detection of all N 
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FIG. 2: (color online). The sensitivity ^mAfl (black solid 
line) for N — 100 calculated with Eq.© as a function of 
\ipin) G A. The sub shot-noise sensitivity is reached for all 
states with non-negligible NOON component (see Eq.lJT]) for 
details). The values of y/rn/S.9sN and y^mA^HL are denoted 
by the upper and lower dashed blue lines, respectively. The 
optimal sensitivity, given by the inverse of the QFI, is drawn 
with the red open circles. 



atoms. It can be demonstrated that if just one particle is 
missed, the sub shot-noise sensitivity is inevitably lost. 

Conclusions. We analyzed different phase estimation 
strategies based on the position measurement of atoms 
released from a double-well trap. We demonstrated that 
the fit to the density gives a sensitivity limited by the 
shot-noise. This bound can be overcome by phase esti- 
mation with correlation functions of the order of at least 
y/N, while the iV-th order correlation function provides 
an optimal estimation strategy. We also showed that the 
measurement of the position of the center of mass gives 
sub shot-noise sensitivity for all states with non-negligible 
NOON component. Both the protocol involving high- 
order correlations and the one based on the center of 
mass position are difficult to implement. The difficulty 
in achieving sub shot-noise sensitivity can be attributed 
to the fact that, after formation of the interference pat- 
tern, the two spatial modes, corresponding to the arms 
of the interferometer, cannot be distinguished and the 
useful information about the correlations between these 
modes is unavailable. To beat the shot-noise limit, one 
has to determine the correlations between the particles. 
These correlations are difficult to extract from the ex- 
perimental data, thus sub shot-noise sensitivity with two 
interfering BECs might prove challenging. 
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